How to make semiconductors ferromagnetic: A first course on spintronics 



m 
o 
o 

(N 

Ph. 
< 



C/2 



o 



> 

O 

cn 

o 



I 

O 

o 



S. Das Sarma, E. H. Hwang, and A. Kaniinski 

Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111 

The rapidly developing field of ferromagnetism in diluted magnetic semiconductors, where a semi- 
conductor host is magnetically doped by transition metal impurities to produce a ferromagnetic 
semiconductor (e.g. Gai-^Mn^^As with a; ~ 1 — 10%), is discussed with the emphasis on elucidating 
the physical mechanisms underlying the magnetic properties. Recent key developments are summa- 
rized with critical discussions of the roles of disorder, localization, band structure, defects, and the 
choice of materials in producing good magnetic quality and high Curie temperature. The correlation 
between magnetic and transport properties is argued to be a crucial ingredient in developing a full 
understanding of the properties of ferromagnetic semiconductors. 



I. INTRODUCTION 

The study of semiconductors (i.e. intrinsically insu- 
lating materials with band gaps in the ^ 1 eV range) 
which are also ferromagnetic has a long, but perhaps 
somewhat intermittent, history. Back in the 1960s and 
early 1970s there was a great deal of interest-'^^ in Eu-based 
chalcogenide materials which exhibited ferromagnetism 
with Curie temperature around 50K or less. These 
were strongly insulating systems with rather poor semi- 
conducting transport properties. In the 1980s there 
were major efforts^ in magnetically doped II- VI ternary 
and quarternary semiconductor alloys, e.g. CdMnSe, 
CdMnTe, PbSnMnTe etc., but most of these materials 
exhibited spin glass or related disordered magnetic be- 
havior and any ferromagnetism invariably had very low 
(< 5K) transition temperature. The idea of a diluted 
magnetic semiconductor, where a small (< 10%) con- 
centration of magnetically active atoms (most often Mn 
although Co and Cr have also been used occasionally) 
is distributed at the cation sites of the host semiconduc- 
tor, is appealing because such a system may have both 
semiconducting and ferromagnetic properties. Making 
such materials has, however, been problematic because 
the magnetic transition metal (i.e. Mn) is not thermody- 
namically stable in the semiconductor host (e.g. GaAs), 
and tends to segregate. 

The need for actively controlling and manipulating 
the carrier spin degree of freedom (in addition to the 
charge degree of freedom) in semiconductors leading 
to novel multifunctional magneto-opto-electronic capa- 
bilities has created the emerging interdisciplinary field 
of spintronics (or spin electronics). A key compo- 
nent of spintronics is the development of new ferromag- 
netic semiconductors. Following the successful devel- 
opment of Gai_a;Mn3:As2i^ and Ini-^^Muj^Aa^ as fer- 
romagnetic semiconductors (with a; « 1 — 10%) using 
careful low-temperature molecular beam epitaxy (MBE) 
technique, intensive worldwide activity has led to claims 
of ferromagnetism (some at room temperatures and 
above) in several magnetically doped semiconductors, 
e.g. GaMnPfL GaMnN,* GeMn,^ GaMnSb.^" It is at 
present unclear whether all these reports of ferromag- 



netism (particularly at room temperatures or above) are 
indeed intrinsic magnetic behavior or are arising from 
clustering and segregation effects associated with vari- 
ous Mn-complexes (which have low solubility) and re- 
lated materials problems. The observed ferromagnetism 
of GaMnAs is, however, well-established and is univer- 
sally believed to be an intrinsic diluted magnetic semi- 
conductor (DMS) phenomenon. 

The Mn dopants in GaMnAs serve the dual roles of 
magnetic impurities providing the local magnetic mo- 
ments and of acceptors producing, in principle, one hole 
per Mn atom. The number density ric of the charge car- 
riers (holes in GaMnAs), however, turns out to be by 
almost an order of magnitude lower than the number den- 
sity Ui of the Mn ions, making the system to be a highly 
compensated doped semiconductor. This high level of 
compensation is believed to be a crucial ingredient in the 
underlying DMS ferromagnetic mechanism. In particu- 
lar, the precise role played by the relative values of the 
local moment density rii and the hole density Uc in giving 
rise to DMS ferromagnetism is currently being debated 
in the literature to the extent that there is no agreement 
even on whether the heavy compensation [ric <^ n{) in 
the system helps or hinders ferromagnetism. An impor- 
tant question in this context is to obtain the functional 
dependence of the ferromagnetic transition temperature 
( "Curie temperature" ) Tc on magnetic moment and hole 
densities rij and ric- In Gai_:rMn2;As the early experi- 
mental resulta^iii indicated a nonmonotonic behavior of 
Tc{x) with Tc first increasing with the Mn concentra- 
tion X, reaching a maximum of about llOK for x « 5%, 
and then decreasing with further increase of Mn con- 
centration above this optimum value. A strong correla- 
tion between the magnetic transition temperature and 
transport properties of GaMnAs was also reported in 
these findingsiii with the system exhibiting relatively 
stronger metallic behavior close to this optimum Mn dop- 
ing level (~ 5%). Recent experimental studiesiSiH of 
Gai-xMn^^As under carefully controlled annealing con- 
ditions seem to suggest that the metastable nature (and 
high defect content) of low temperature MBE grown 
GaMnAs may be playing an important role in determin- 
ing magnetic properties - in particular, careful annealing 
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could, under suitable conditions, lead to an enhancement 
of Tj, with increasing Mn content leading to an eventual 
saturation of the transition temperature in the x — 5—8% 
range without the nonmonotonic Tc(x) behavior observed 
earlier. These annealing studiesiSii^ confirm the close 
connection between carrier transport properties and fer- 
romagnetism in DMS materials showing that typically 
Tc values are the highest in GaMnAs samples with nomi- 
nally highest conductivities (and presumably the highest 
hole densities). It has therefore been suggested that, in 
agreement with the mean-field theory (MFT) for DMS 
ferromagnetism originally derived in Refs. |]3l5l in a dif- 
ferent context and recently rediscovered in Refs. H^IT^ 

1 /3 

Tc in GaMnAs grows monotonically as Tc oc nj when 
the carrier density increases for a fixed Mn concentra- 
tion. This simple Weiss MFTili^JMLi^ then implies 
a convenient technique for enhancing the ferromagnetic 
transition temperature in GaMnAs, namely, to increase 
the hole density as much as possible. The general validity 
of this Weiss MFT-based assertion remains questionable, 
and it has been arguedi2iS£ that this prediction holds 
only in perturbative regime of small J and small Uc- 

The GaMnAs annealing studies bring up important 
and interesting questions regarding the nature and role of 
metastable defects in the system that are presumably be- 
ing healed during annealing leading to Tc enhancement. 
Two different kinds of GaMnAs lattice (atomistic) de- 
fects have been discussed in this context. First, there are 
As antisite (AsQa) defects. As atoms sitting at Ga sites, 
that are invariably present in the low temperature MBE 
growth required for DMS GaMnAs. These As antisite de- 
fects act as double donors with each As defect capturing 
two holes in GaMnAs. The second type of defects that 
has been much discussed recently, is the Mn interstitial 
defect, where the "defective" Mn atom sits at an intersti- 
tial site in the GaAs lattice instead of occupying a sub- 
stitutional cation site replacing a Ga atom. The Mn in- 
terstitials are thought to be particularly deleterious in af- 
fecting the magnetic properties of GaMnAs because they 
have two adverse effects: first, they act as donors (similar 
to As antisite defects) capturing holes; second, the Mn 
interstitials typically have strong local direct antiferro- 
magnetic coupling to any neighboring substitutional Mn 
atoms, thereby suppressing the net effective Mn magnetic 
moment available for DMS ferromagnetism. The quan- 
titative abundance of As antisites and Mn interstitials 
in GaMnAs is difficult to estimate because there is really 
no direct experimental technique to observe these defects. 
It is, however, widely believed that the low-temperature 
MBE technique, crucial for the homogeneous single-phase 
growth of GaMnAs, invariably generates both of these 
defects with their relative importance varying from sam- 
ple to sample depending on the uncontrollable details of 
growth and processing. In addition to the atomistic As 
antisite and interstitial defects in GaMnAs, there are in- 
variably many other unintentional (as well as unknown 
and uncontrolled) defects and impurities (both charged 
and neutral) in GaMnAs because low temperature MBE, 



in general, produces materials of rather poor quality. 
This poor quality of GaMnAs is reflected in the trans- 
port properties as even the best GaMnAs samples seem 
to have very low carrier mobilities. The measured hole 
mobilities in the so-called "metallic" GaMnAs systems 
are only about 1 — 50 cm'^/V-s corresponding to a typ- 
ical carrier mean free path of 0.1 — 1 nm. Such small 
values of mean free paths would indicate a strong ten- 
dency toward carrier localization in GaMnAs, and indeed 
there are reports of reentrant metal-insulator transitions 
in GaMnAs as a function of the Mn content x. Specif- 
ically, it has been reportedii that Gai-^iMn^^As is an 
insulator for small (x < 3%) and large [x > 7%) values 
of X whereas it is metallic (albeit with very small mean 
free paths) in the intermediate range of x (~ 5%) where 
the ferromagnetic transition temperature is maximized. 
Other experimental studies (often using careful sample 
annealing procedures) do not find such reentrant metal- 
insulator transition behavior. ^^■^■^ 

It is important to emphasize that the existence of DMS 
ferromagnetism seems to be independent of the system 
being metallic or insulating. For GaMnAs, both metal- 
lic and insulating samples are ferromagnetic with the 
transition temperature typically being higher for metallic 
systems although this may not necessarily be a generic 
behavior. Many other systems exhibiting DMS ferro- 
magnetism, e.g., InMnAs, GeMn, GaMnSb, GaMnN, are 
however strongly insulating, and as such, metallicity (i.e. 
itinerant nature of the carriers), is not a necessary pre- 
condition for the existence of ferromagnetism in DMS 
systems. Additionally, in Gai-^^Mn^jAs samples exhibit- 
ing reentrant metal-insulator transitions the ferromag- 
netic behavior appears to be completely continuous as a 
function of x with the only observable effect being a vari- 
ation in Tc, which is expected since different GaMnAs 
samples, even with nominally the same value of x, show 
wide variations in Tc values. The real DMS carriers me- 
diating the ferromagnetic interaction between the local 
moments are generically likely to be far from free holes 
in the valence band of the host semiconductor material. 
They are, in all likelihood, extended or bound carriers 
in the impurity (or the defect) band which forms in the 
presence of the Mn dopants. There is a great deal of 
direct experimental support for the relevance of this im- 
purity band picture for DMS ferromagnetism, and first 
principles band structure calculations also confirm the 
impurity band nature of the carriers active in DMS fer- 
romagnetism. Both the impurity band nature of the car- 
rier system and the presence of strong charge and spin 
disorder in the system (due to the random Mn dopants 
as well as other defects, disorder, and impurities invari- 
ably present in the DMS material) make carrier local- 
ization a relevant issue in DMS ferromagnetism. We be- 
lieve that understanding the correlation between trans- 
port and magnetic properties is a key ingredient in DMS 
physics. 

In the next section of this review we critically discuss 
some of the proposed theoretical models (and their va- 
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lidity or applicability) for DMS ferromagnetism empha- 
sizing the connection between magnetization (as well as 
the Curie temperature) and transport properties. We 
conclude in section III with a brief discussion of (many) 
open questions and some possible future directions in the 
subject. 

This article is is more of a theoretical perspective; it is 
not meant to be an exhaustive review of the vast litera- 
ture on ferromagnetic semiconductors. We only discuss 
some of the key developments in the subject in order 
to highlight the physical mechanisms and the important 
system parameters controlling ferromagnetism in diluted 
magnetic semiconductors. 



II. THEORIES 

The fact that the DMS arc ferromagnetic independent 
of their weakly metallic or strongly insulating nature im- 
plies a robust character for the underlying mechanism 
leading to the long range magnetic order in these sys- 
tems. Clearly the ferromagnetic mechanism should not 
depend crucially on the carrier system being "free" va- 
lence band holes since the strongly insulating DMS sys- 
tems do not have any free holes. The currently accepted 
picture for DMS ferromagnetism is that it is the local an- 
tiferromagnetic coupling between the carriers (i.e., holes 
in GaMnAs) and the Mn magnetic moments that leads 
to long range ferromagnetic ordering of Mn local mo- 
ments. The carrier system also becomes spin-polarized 
in the process with the carrier magnetic moment directed 
against the Mn magnetic ordering by virtue of the anti- 
ferromagnetic hole-Mn coupling, but the total magnetic 
moment of the spin polarized carriers is extremely small 
since ric <C ni and |S| > |s| where S and s are respectively 
the Mn and the hole spin. The relevant DMS effective 
magnetic Hamiltonian can be written as 



Hm^ y"dVj(r)S(r).s(r), 



(1) 



where S(r) and s(r) are respectively the Mn and hole spin 
densities. The coupling J(r) between Mn local moments 
and hole spins can, in principle, be ferromagnetic (J < 0) 
or antiferromagnetic ( J > 0), but the effective interaction 
between the Mn local moments mediated by the holes 
(through Hm) is always ferromagnetic. The magnitude 
of J must come from a first principles band structure 
calculation or from experiments. 

Since the Mn moments are localized at specific lattice 
sites, it is sensible to write Eq. ^ in the more conven- 
tional form 



M 



d3r5(r-R,)S,-s(r) , 



(2) 



where the sum over i goes over the Mn sites in the GaAs 
lattice, and we have assumed the exchange coupling J to 
be a constant. The JS^-s local exchange coupling defined 



in Eq. Q is sometimes referred to as the Zener model 
(or the s — d coupling although in the context GaMnAs, 
where the carriers are holes, the coupling is more like an 
p — d exchange coupling between Mn d levels and the p 
valence band). 

The magnetic Hamiltonian defined in Eqs. Sind 
(|^ is, of course, one of the most well-known interac- 
tion Hamiltonians in all of condensed matter physics. 
Depending on the relative magnitudes of hole and Mn 
moment densities ric and rij, and the coupling strength 
J, Hi\i describes the Kondo model or the Kondo lattice 
model or s — / (or s — d) Zener model or the double ex- 
change model or the RKKY model or (in the presence 
of strong disorder) the (RKKY) spin glass model. Most 
of the earlier studies of this exchange Hamiltonian have 
concentrated on the low impurity-density limit, rii <C ric, 
where the local moments are a perturbation on the car- 
rier Fermi surface. The DMS systems provide an inter- 
esting novel regime for studying Hm, namely the regime 
of low carrier density nj 3> tt-c • There has not been much 
earlier work focused on the low carrier-density limit of 
the Kondo lattice Hamiltonian or the RKKY spin glass 
Hamiltonian, which is the appropriate parameter regime 
for GaMnAs. We emphasize that Kondo physics plays no 
essential role in DMS systems, but because of the strong 
inherent disorder in these systems spin glass physics may 
very well be playing a role at lower values of ni/ric. 

The simplest way to understand DMS ferromagnetism 
on a qualitative level is to neglect all disorder effects, and 
assume that the system can be thought of as a collection 
of local moments of density rij interacting with itiner- 
ant holes of density ric- In the weak coupling regime 
(J <C t, Ep where t and Ep are the band width and the 
Fermi energy of the carrier system) it is then possible 
to "eliminate" the carrier degrees of freedom by map- 
ping the problem into the corresponding RKKY prob- 
lem of Mn local moments interacting indirectly via the 
holes through Hm- This carrier-polarization-mediated 
RKKY exchange interaction is essentially a second order 
perturbation theoretic calculation in J, leading to the 
following interaction Hamiltonian between the Mn spins 
(assuming, for the sake of simplicity, a single parabolic 
hole band) 



Hi — ^ JijSi-Sj, 



with 



sm{2kFR) - {2kFR) cos{2kFR) 



{2kFRr 



(3) 



(4) 



1/3 J. ^ 

where kp oz ric oc Ep" is the Fermi wave vector for the 
holes and R = |Ry| = |Ri — Rj | is the distance between 
the Mn impurities at i and j lattice sites. This oscillatory 
RKKY indirect exchange magnetic coupling between lo- 
cal moments placed in an electron gas arises from the spin 
polarization induced in the electron system by the local 
moment, and is a perturbative effect of the local moments 



,1/2 
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on the carrier system (valid in the hmit of small J). In 
the DMS systems an approximate average distance be- 

~ 1/3 

tween the Mn local moments is R , making typical 

values of fci?i? '--^ (nc/"-i)^^^- For ferromagnetism to occur 
(remembering that the real GaMnAs system is intrinsi- 
cally strongly substitutionally disordered) one must avoid 
the strong frustration effects associated with random sign 
changes in the magnetic interaction of Eq. ((2Jl arising 
from the oscillatory terms in Eq. which can only be 
guaranteed when kpR <C 1. This then necessitates the 
condition ric/ni 1, making strong compensation, i.e., 
ric <^ Hi, a necessary condition for ferromagnetism in the 
system. 

We note that a Weiss MET calculation can easily be 
carried out for the coupled Mn-hole magnetic problem 
defined by the Zener-RKKY model of Eqs. CEll-O. Eor a 
simple single-band degenerate hole model, the result for 
in the Weiss MET was obtained a very long time 

rpMF 5(5+1) 2|7|2n/E^ \ fr^\ 

= ^ s \ J\ D[EF)ni, (5) 

where S is the local moment spin and s is the carrier 
spin. The carrier density of states D{Ep) in Eq. jSJ, 
in the degenerate one band free carrier system, is given 

by D{Ep) ex cx 71^^. Thus increases mono- 

tonically with the exchange coupling, the local moment 
density, and the carrier density as 

T^" oc J^nl'^n,, (6) 

where Uc and n-i are treated as independent density pa- 
rameters. Based on this simple Weiss MET (and its 
various extensions incorporating spin wave corrections 
and/or explicit spin-orbit coupling effects in the GaAs va- 
lence band), quantitative predictions have been made^® 
about the dependence of the ferromagnetic transition 
temperature in various DMS systems on the system pa- 
rameters (carrier and local moment densities, effective 
exchange coupling strength, etc.). In addition, we have 
recently shown^^ that the MET for the model defined by 
Eqs. ^ and |(2J) can be extended to the nondegenerate 
carrier system more appropriate for the localized holes 
in the strongly insulating DMS materials. In contrast to 
the degenerate results of Eq. ||HJl, which is appropriate 
for the metallic ("extended") carrier system (e.g. free 
valence band holes), the nondegenerate MET provides a 
transition temperature dependence given by 

^\J\{n,n,Yl^. (7) 

Thus, the dependence of on the system parameters 
changes drastically in the localized nondegenerate limit. 
More details can be found in Ref. |^ 

The simple Weiss mean- field theory (Eqs. (O and ©) 
based on the model of degenerate valence band holes in- 
teracting with the Mn local moments (and its obvious 
extensions) has been used extensively in the literature. 



The question, therefore, arises about its validity, or more 
generally, about its quantitative reliability. It is argued 
that the molecular MET should apply well in DMS sys- 
tems because of the large length scales (e.g. fc^^) involved 
in the problem. We believe that there are good reasons 
to question the validity of the Weiss MET as applied to 
DMS ferromagnetism. Eirst, strong disorder in the sys- 
tem invariably brings in some effects of frustration which 
should become more significant at larger values of car- 
rier density, making the simple MET defined by Eqs. Q 
and (jSJ increasingly incorrect at higher carrier densities 
(even for the metallic system). In fact, Tc should vanish 
in the disordered system as Uc Ui. Second, MET in 
this problem being equivalent to a second order pertur- 
bation theory in the exchange coupling J (i.e. Eqs. QJ- 
©), the validity of Weiss mean-field predictions becomes 
highly questionable for large values of J where carriers 
may bind to individual local moment sites to lower their 
energy forming an impurity band. There is significant ex- 
perimental evidence for impurity band formation in DMS 
materials, which is not included in the MET. Therefore 
mean-field predictions will fail at larger values of J where 
DMS magnetic properties should be dominated by impu- 
rity band physics not caught by Eqs. Q and ©. We 
have recently developed a dynamical mean-field theory 
(DMET) picturai2iSS of DMS ferromagnetism where im- 
purity band physics is explicitly included. Third, in the 
limit of rtj ^ Uc, where DMS ferromagnetism prevails, 
disorder effects are extremely strong (as is obvious from 
the insulating nature of most DMS materials and from 
the extremely low mobilities of the so-called "metallic" 
DMS materials which may very well also be insulating 
at T = 0), and therefore the degenerate "metallic" Weiss 
MET approach based on the carriers being the free ex- 
tended valence band holes of the host semiconductor may 
be fundamentally flawed. A more appropriate picture in 
this disordered localized regime may be the percolation 
picture of DMS ferromagnetism recently developed by 
us^ where each hole is bound in a cluster of magnetic 
impurities forming a bound magnetic polaron. With de- 
creasing temperature, these bound magnetic polarons co- 
alesce eventually leading to a magnetic percolation tran- 
sition at Tc where a magnetic percolation cluster spans 
the whole system. This polaron percolation picture of 
DMS ferromagnetism is strongly supported^^ by direct 
numerical simulations and by experimental studies show- 
ing strongly concave magnetization curves in many DMS 
systems which follow naturally from our theoryi2iiS^ 

Below we provide some details on the dynamical mean- 
field and polaron percolation theories of DMS ferromag- 
netism. 



A. Dynamical mean-field theory 

In this section we briefly describe the "dynamical 
mean- field theory" (DMET). Due to size limitations, we 
are not able to give detailed description of this formalism 
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here. A reader i nterested in more details may find them 
in Ref. 11912(1211 

We model the Gai-^^Mna; As system as a lattice of sites, 
which are randomly nonmagnetic (with probability 1 — x) 
or magnetic (with probability x), where x now indicates 
the relative concentration {i.e. per Ga site) of active Mn 
local moments in Gai^^jMn^jAs. The DMFT approxima- 
tion amounts to assuming that the self energy is local or 
momentum independent, S(p, zw^) S(ia;„), and then 
all of the relevant physics may be determined from the 
local (momentum-integrated) Green function defined by 

/d^p 1 
7T^' — [ n — — ^' 
[2ny^ lUn + M - e(P) - ^a(lUJn) 

where we have normalized the momentum integral to the 
volume of the unit cell Oq, and is the chemical poten- 
tial. Gioc is in general a matrix in spin and band indices 
and depends on whether one is considering a magnetic 
or non-magnetic site. Since Gioc is a local function, it is 
the solution of a local problem specified by a mean-field 
function go, which is related to the partition function 
Zioc = / dS cxp(— 5100) with action 

Sloe = J29Lpir^r')ci{r)c^,{r') 

+ j5](S.^„^)ct(r)c^(r'), (9) 

a/3 

on a magnetic site and 

5ioc = EfoaM^ - ^')ci(r)c^(r'), (10) 

on a non-magnetic site. Here c^{t) (c]j(t)) is the destruc- 
tion (creation) operator of a fermion in the spin state a 
and at time r. go{T — T') plays the role of the Weiss mean 
field (bare Green function for the local effective action 
5ioc) and is a function of time. The magnetic-site mean- 
field function can be written as 5oa/3 = ao + aim.- (Tap 
with m the magnetization direction and ai vanishing in 
the paramagnetic state (T > Tc). Since the spin axis 
is chosen parallel to m, becomes a diagonal matrix 
with components parallel {g^-^ = + ai) and antipar- 
allel [g^^ = ao — oi) to m. It is specified by the condi- 
tion that the local Green function computed from Zioc, 
namely JlnZioc/^ffo = (ffg ~ ^)~^ is identical to the lo- 
cal Green function computed by performing the momen- 
tum integral using the same self energy. 

Within this approximation we calculate the normalized 
magnetization of the local moments with the semicircle 
density of states, D{uj) = \f^\?^^Zp- jl-Kt, where t is a 
bandwidth of the system. The normalized magnetization 
M(r) is given by M(T) = / (dS) • S exp(-5ioc)/^ioc. As 
the temperature is increased, the spins disorder and even- 
tually the magnetic transition temperature is reached. 
Above this temperature, 50 is spin- independent. By lin- 
earizing the equation in the magnetic part of go with 




FIG. 1: Calculated ferromagnetic transition temperature 
as a function of carrier carrier density n^jx for x — 0.05 
and various exchange couplings J. Inset shows the majority 
spin density of states aX T = and for x = 0.05 and J = 
0.5t, l.Qt. 

respect to ai we may obtain the ferromagnetic transition 
temperature Tn. The details on the calculation of Tc are 
given m Refs. Il9l2nt In Fig. □ we show the calculated 
Tc as a function of carrier density for x — 0.05 and vari- 
ous couplings J, and the corresponding density of states 
are given in the inset. When J < t the impurity band 
formation dose not occur (see the inset for J = 0.5i) 
and Tc are saturated as the density increases. However, 
as J > t the impurity spin bands are well established 
(see the inset for J = 2.0t). In this limit the Tc has a 
maximum value around the half filling of the impurity 
band [uc ~ O.bx) and Tc vanishes as the carriers fill the 
band {uc = x) because no low energy hopping processes 
are allowed in a ferromagnetic state for a filled impurity 
band. Within DMFT the calculated Tc shows the non- 
momotonic behavior for all J and the Weiss mean field 
prediction Tc oc n}/^ can be applied only for small Uc- 

In Fig. 121 we show the normalized magnetization of the 
local moments as a function of temperature for differ- 
ent values of J =1.0, 1.5, 2.Qt and x — 0.05, and for 
various hole densities, ric/rti. For the small coupling con- 
stant J — t the impurity band is not formed, but for 
J =1.5, 2.0t we have a spin polarized impurity band. 
For relatively high density {uc/rii = 0.4) the calculated 
magnetization looks similar to the Weiss mean-field re- 
sults. But for low density (rtc/«i = 0.04) we have a linear 
M{T) in the intermediate temperature range. Near the 
critical temperature Tc the critical behavior of the mag- 
netization for aU density is given by M{T) oc (Tc — T)^/^. 
For different exchange couplings we have similar results 
(i. e., linear behavior at low densities and intermediate 
temperature ranges). 
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FIG. 2: The normalized DMFT impurity magnetization as a function of temperature for (a) J = 1.5t, (b) J = l.Qt, and (c) 
J = 2.0t, for a; = 0.05 and for ric/rii =0.4, 0.2, 0.1, 0.04 (from the top). The dashed line in (a) represents the magnetization 
calculated for the simple Weiss mean field theory for the local moment spin S = 5/2. 



In particular, it should be possible to obtain the M (T) 
behavior for the localized carrier case also from DMFT by 
incorporating impurity band localization in the DMFT 
formalism. Our current theory does not include localiza- 
tion, and the impurity band (or valence band) carriers in 
our DMFT calculations are all delocalized metallic car- 
riers. First principles band theory calculations indicate 
that the actual exchange coupling in Gai-ajMn^^As may 
be close to critical Jc, and as such impurity band physics 
may be quite important for understanding DMS magne- 
tization. To incorporate physics of localization one needs 
to include disorder effects (invariably present in real DMS 
systems) in the model. All our DMFT calculations are 
done in the virtual crystal approximation where effective 
field is averaged appropriately leaving out random dis- 
order effects explicitly. In the next section we explicitly 
incorporate disorder in the theory by developing a per- 
colation theory approach to DMS magnetization for the 
strongly localized insulating systems. 



B. Percolation theory 

Our percolation theory for DMSSS applies strictly in 
the regime of strongly localized holes where the dynam- 
ical mean-field theory for delocalized carriers described 
above has little validity. These two theories, mean-field 
theory and percolation theory, are therefore complemen- 
tary. 

The percolation theory assumes the same model of 
carrier-mediated ferromagnetism, but now the carriers 
are pinned down with the localization radius as- The 
disorder, averaged out in the mean-field theory, plays a 
key role in the carrier localization. The magnetic impuri- 
ties are assumed to be completely randomly distributed 
in the host semiconductor lattice in contrast to the mean- 
field case where the carrier states are free and the disorder 
is neglected. 

As we have demonstrated in Ref. [23 . the problem of 
ferromagnetic transition in a system of bound magnetic 



polarons can be mapped onto the problem of overlap- 
ping spheres well-known in the percolation theory. The 
latter problem studies spheres of the same radius r ran- 
domly placed in space (three-dimensional in our case) 
with some concentration n. Overlapping spheres make 
"clusters;" as the sphere radius r becomes larger, more 
and more spheres join into clusters, the clusters coalesce, 
and finally, at some critical value of the sphere radius, an 
infinite cluster spanning the whole sample appears. This 



problem has only one dimensionless parameter. 



r n. 



and 



therefore can be easily studied by means of Monte-Carlo 
simulations. 

Each sphere of the overlapping spheres problem corre- 
sponds to a bound magnetic polaron, which is a complex 
formed by one localized hole and many magnetic impuri- 
ties with their spins polarized by the exchange interaction 
with the hole spin. The concentration n of spheres coin- 
cides with the concentration ric of localized holes. The 
expression for the effective polaron radius in terms of the 
physical parameters of the system under consideration is 
not trivial and has been found in our earlier workiSS 



r n 



0.86 + {alric) " In 



(11) 



Here 0.64 « 0.86"^ is the critical value of the parame- 
ter r^n at which the infinite cluster appears, and Tc is 
the Curie temperature of the ferromagnetic system under 
consideration, derived in Ref. i22, 



Tr. 



as 



exp 



0.86 

(alnc)' 



(12) 

The limit of applicability of Eq. (|12|l is determined by 
the condition a^ric <C 1. The dependence of Tc on a^n^ 
is shown in Fig. (21 

Since Uc <C ni, one polaron includes many magnetic 
impurities, and the total magnetization of the sample is 
that of impurities: 



M{T) 
Mo 



= V 0.86+ (airic 



T 



(13) 
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FIG. 3: Curie temperature Tc as a function of the dimension- 
less parameter a%nc. At a%nc < 1, Tc is given by Eq. 11211 . At 
a%nc > 1, Eq. I|12^ (being beyond limits of its applicability ) 
predicts decline of Tc (dotted line); in reality, Tc grows mono- 
tonically with a%nc (solid line), though its exact behavior is 
unknown. 

where the universal function V{y) is the infinite cluster's 
volume in the model of overlapping spheres; it depends 
only on the product y of the spheres' diameter and the 
cubic root of their concentration. One can see that the 
shape of the magnetization curve is determined by only 
one dimensionless parameter a^ric, while the expression 
(I12|l for Tc is more complicated and depends on all pa- 
rameters of the model. Fig. 01 shows the temperature 
dependence of the magnetization at two values of a^ric; 
the curve is more concave at smaller values of this param- 
eter. The concave shape of the curve is consistent with 
the experimental magnetization data in Gei-ajMn^;^ and 
low-Tc ni-V samples^' where our polaron percolation 
picture applies better due to stronger carrier localiza- 
tion associated with lower values of a^ric- Our magne- 
tizat ion re sults also agree with the numerical results of 
Ref. I23I25I Very recently, a direct numerical Hartree- 
Fock calculation has convincingly verified the magnetic 
polaron percolation picture for the localized regimei^^ 



III. CONCLUSION 

When carriers are present with a carrier density ric 
in a semiconductor containing a certain magnetic im- 
purity density rii of local magnetic moments, a Weiss 
mean field theory, developed years ago,^^'^^ within the 
Zener-RKKY effective magnetic Hamiltonian describing 
the local exchange coupling between local moments and 
carrier spin, would always predict ferromagnetism with a 
mean field transition temperature varying linearly 

1 /3 

in Hi and Uc ■ This Weiss mean field theory, appro- 
priately modified by the complicated spin-orbit coupled 
multiband valence band structure of GaAs, has been ex- 
tensively used in making quantitative predictions about 
the magnetic properties of GaMnAs. This Weiss mean 




FIG. 4: The solid lines show the relative magnetization of 
the magnetic impurities [Eq. for a%nc — 10~^ (curve 

1) and 10~^ (curve 2). The dashed lines show the relative 
magnetization of localized holes, whose contribution to the 
total sample magnetization is small. 



field theory, while being a powerful qualitative tool, is 
obviously too simplistic to be taken seriously for quanti- 
tative analyses of DMS properties. In real DMS systems 
spatial fluctuations (i.e. disorder associated with random 
positions of the magnetic dopants) and thermal fluctua- 
tions of magnetic moments as well as impurity band and 
discrete lattice effects (all neglected in the Weiss mean 
field theory) are likely to be of real importance even if the 
magnetic coupling can be assumed to be a simple local 
exchange coupling between local impurity moments and 
carrier spins. In addition, the mean field theory averages 
out the oscillatory exchange coupling between the impu- 
rity moments, which must be important at higher carrier 
densities, particularly in the presence of strong spatial 
disorder. For example, the mean field theory incorrectly 
predicts the ground state to be ferromagnetic for all val- 
ues of ric and rii, whereas in reality of course the system 
is a spin glass for larger values of ric. The mean field the- 
ory is qualitatively valid only for very small values of the 
carrier-local-moment exchange coupling J where impu- 
rity band effects are weak, and applies only in the limit 
of Tie <^ Hi where disorder and oscillatory aspects of ex- 
change can be neglected. In addition, the standard mean 
field theory is valid only for degenerate metallic systems 
although a nondegenerate localized version of the mean 
field theory has recently been developed. 

There have been recent theoretical 
development3i2i22i2iiSS going beyond the simple Weiss 
mean field theory where impurity band effects (both 
for metalli(}i2i2£ and insulating-^ DMS systems) are 
explicitly included in the dynamical mean fieldiSiSfliSi 
and magnetic polaron percolation theoriesSiiSS. The 
DMFT theory predicts, in contradiction with the Weiss 
mean field theory, that the ferromagnetic transition 
temperature Tc is not monotonic in the carrier density 
(and the coupling strength J), and may vanish for 
larger values of Uc- While the definitive quantitative 
theory for DMS ferromagnetism, taking into account 
complete band structure complications of the system 
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along with disorder, fluctuation, and nonperturbative 
strong-coupling (i.e. J not necessarily small) effects in a 
consistent manner is still lacking, it is increasingly clear 
that disorder and impurity band aspects of the real DMS 
systems are extremely important in understanding DMS 
magnetic properties, and theoriesi^iSLSiiS^ that include 
these effects should be of qualitative and quantitative 
significance in explaining DMS magnetic properties. 

The observed correlation between transport and mag- 
netic properties in DMS systems, e.g. the more metal- 
lic systems typically having more convex magnetiza- 
tion curves as a function of temperature whereas more 
insulating systems having extreme non-mean-field like 
concave magnetization versus temperature curves (with 
considerable missing magnetic moments), is perhaps a 
strong indicator of the underlying mechanism of DMS 
ferromagnetism.^^ This correlation indicates a carrier- 
mediated magnetic interaction origin of DMS ferro- 
magnetism. In particular, two very recent theoretical 
publication32iiS& have discussed in considerable details 
this correlation between metallicity and DMS ferromag- 
netism. In Ref. |^ it has been argued that the appro- 
priate picture for DMS ferromagnetism in strongly disor- 
dered insulating materials is the bound magnetic polaron 
percolation picture whereas at higher densities (and/or 
lower disorder), where the system is metallic, the dynam- 



ical mean-field picture (which becomes equivalent to the 
Weiss mean field picture at low values of exchange cou- 
pling insufhcient to cause the impurity band formation) is 
more appropriate. In Ref. 26 a detailed numerical simu- 
lation of the disordered DMS model has been carried out 
qualitatively validating this picture. The question there- 
fore arises regarding the nature of the DMS ground state: 
Is it localized and insulating or extended and metallic or 
is there a transition between the two as a function of 
carrier density? If there is a true metal-insulator transi- 
tion, then the percolation picture holds in the insulating 
phase and the dynamical mean field picture holds in the 
metallic phase. Although there have been experimental 
claims of metal-insulator transitions in GaMnAs systems 
as a function of Mn content, the current experimental 
situation is sufficiently murky that is possible that all 
GaMnAs systems are in fact insulating at T = for all 
values of the Mn content. It would therefore be desir- 
able to develop a dynamical mean field theory (or even 
just the Weiss mean field theory) of DMS ferromagnetism 
explicitly incorporating disorder effects so that the mag- 
netic properties can be studied continuously through the 
metal-insulator transition (or crossover, as the case may 
be). 

This work is supported by US-ONR and DARPA. 
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